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A clear demonstration of topological superconductivity (TS) and Majorana zero modes remains 
one of the major pending goal in the field of topological materials. One common strategy to generate 
TS is through the coupling of an s-wave superconductor to a helical half-metallic system. Numerous 
proposals for the latter have been put forward in the literature, most of them based on semiconduc¬ 
tors or topological insulators with strong spin-orbit coupling. Here we demonstrate an alternative 
approach for the creation of TS in graphene/superconductor junctions without the need of spin- 
orbit coupling. Our prediction stems from the helicity of graphene’s zero Landau level edge states 
in the presence of interactions, and on the possibility, experimentally demonstrated, to tune their 
magnetic properties with in-plane magnetic fields. We show how canted antiferromagnetic ordering 
in the graphene bulk close to neutrality induces TS along the junction, and gives rise to isolated, 
topologically protected Majorana bound states at either end. We also discuss possible strategies to 
detect their presence in graphene Josephson junctions through Fraunhofer pattern anomalies and 
Andreev spectroscopy. The latter in particular exhibits strong unambiguous signatures of the pres¬ 
ence of the Majorana states in the form of universal zero bias anomalies. Remarkable progress has 
recently been reported in the fabrication of the proposed type of junctions, which offers a promising 
outlook for Majorana physics in graphene systems. 


I. INTRODUCTION 

The realisation of topological superconductivity (TS), 
a novel electronic phase characterised by Majorana 
excitations, has become a major goal in modern 
condensed matter research. Despite promising ex¬ 
perimental progress pUT3] on a number of appealing 
implementations pMUTT] . a conclusive proof of TS re¬ 
mains an open challenge. We here report on a 
new approach to obtain TS and Majorana states in 
graphene/superconductor junctions. Key to our pro¬ 
posal is the interaction-induced magnetic ordering of 
graphene’s zero Landau level (ZLL). Coupling this unique 
state to a conventional superconductor gives rise to novel 
edge states whose properties depend on the type of mag¬ 
netic order. In particular, the canted antiferromag¬ 
netic phase is a natural host for Majorana bound states. 
Our proposal combines effects that were recently demon¬ 
strated experimentally (tunable spin ordering of the ZLL 
na cog and ballistic [20j [21] graphene/superconductor 
junctions of high-transparency [21] operating in the 
Quantum Hall regime [20]), and is thus ready to be 
tested. 

While intrinsic TS is rare, it can be synthesised ef¬ 
fectively through the coupling of a conventional s-wave 
superconductor (SC) and tailored electronic gases with 
spin-momentum locking. Using this recipe, it has been 
predicted that Majorana excitations should emerge when 
one induces superconductivity onto topological insulators 
[T4] or semiconductors with strong spin-orbit coupling 
p~5j . Particularly attractive are implementations of one¬ 


dimensional TS using either semiconducting nanowires 
mnzi or edge states in two-dimensional Quantum Spin 
Hall (QSH) insulators, since the main ingredients are al¬ 
ready in hand. These ideas have spurred a great deal of 
experimental activity mm- Despite this progress, how¬ 
ever, an unambiguous demonstration of TS is, arguably, 
still missing. Important limitations of these systems in¬ 
clude disorder, bulk leakage, or imperfect proximity effect 
(the so-called soft gap problem). Thus, it is worthwhile 
to explore alternative materials. 

One particularly interesting option is graphene [22], 
which exhibits very large mobilities even in ambient con¬ 
ditions, and where a ballistic proximity effect has been 
recently demonstrated [21 . Graphene was the first ma¬ 
terial where a topological insulating phase was proposed 
[ 23] in the presence of a finite intrinsic spin-orbit cou¬ 
pling. Kane and Mele showed that graphene then be¬ 
comes gapped around neutrality, and a single helical edge 
mode with spin locked to propagation direction develops 
at each edge. In such QSH regime, gapping the edge 
states through proximity to a conventional superconduc¬ 
tor gives rise to a one-dimensional TS along the interface 
m- Graphene’s negligible spin-orbit coupling, however, 
has proved to be a fundamental roadblock in this pro¬ 
gramme. 

In this work we present a simple mechanism to re¬ 
alise the above situation in graphene without recourse 
to spin-orbit coupling. We consider a graphene ribbon 
in the Quantum Hall (QH) regime, in which, unlike in 
the QSH case, time-reversal symmetry is broken by a 
strong magnetic flux (Appendix [A]). In contrast to con- 
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FIG. 1. Sketch, bandstructure and phase diagram of a 350 nm-wide graphene sample (Fermi velocity vf = 10 6 m/s) in the 
Quantum Hall regime (out of plane field B z — 1 T) in various configurations. The chemical potential (dotted line) is tuned 
within the gap Azll ~ 20meV of the zero Landau level, which has ferromagnetic (a-c), antiferromagnetic (d-f) or canted 
antiferromagnetic ordering (g-i) due to electronic interactions. The bandstructures under each sketch correspond to an infinite 
graphene ribbon surrounded by vacuum (vac/Gr/vac panels) or coupled to a superconductor of gap Asc (chosen large ~ 7 meV 
for visibility) along the top edge as in the sketches (vac/Gr/SC panels, Nambu bands). Bands in red correspond to eigenstates 
localised at the top edge of the ribbon. The zero-energy local density of states is shown in blue in each sketch, (j) Low energy 
bandstructure of states along a graphene/superconductor interface folded onto the T point, in the gapless (left panel, zoom of 
panel f), non-trivially gapped (middle panel, canted order) and trivially gapped (right panel, with strong intervalley scattering) 
phases, (k) Phase diagram of said interface, computed from its low energy effective Hamiltonian (see text), as a function of 
magnetic angle 0 and intervalley coupling w. The three possible phases are shown, bounded by threshold values wo. Wins and 
0i ns , see main text. 


ventional two-dimensional electron gases, graphene de¬ 
velops a zero-Landau level at the Dirac point, which has 
been shown in pristine samples to become split due to 
electronic interactions ]24f[29j . Experimental evidence 
m points towards spontaneous antiferromagnetic order¬ 
ing [33, 31], although other broken symmetries have been 
discussed [T8|. In this work we consider all possible mag¬ 
netic orders. Fig. [l] summarises the different possibili¬ 
ties, ranging from ferromagnetic (F) to antiferromagnetic 
(AF) ordering [32] . including canted AF which may be 
controlled by an external in plane Zeeman field as argued 
in Ref. |19] The different orders are parameterised by the 
angle 0 between the spin orientation of the ZLL in the 
two graphene sublattices, so that 6 = 0 for F and 0 = tt 
for AF. 

While 0 is considered in our model as an externally 
tuneable parameter as in Ref. IH we have checked that 
a mean-field calculation in a honeycomb Hubbard model 
under an in-plane Zeeman field (Appendix [A]) yields the 
same bulk and edge phenomenology presented in this 
work m • Corrections beyond mean field and the Hub¬ 
bard model have been explored theoretically in the past, 
and have predicted the formation of a Luttinger liquid do¬ 
main on an infinite vacuum edge [33] . The corresponding 
excitation density resembles the non-interacting edge for 
F order, rather than the AF case. The interacting prob¬ 
lem at a highly transparent superconducting contact re¬ 
mains an open problem. We conjecture that, given their 
robust topological origin, the Majorana phenomena de¬ 


scribed here at a mean field level would survive in the 
Luttinger regime, at least within a limited range of pa¬ 
rameters, and with power-law corrections to the trans¬ 
port results. These issues, however, remain beyond the 
scope of this work. 


II. TOPOLOGICAL SUPERCONDUCTIVITY IN 
QUANTUM HALL GRAPHENE 

An infinite ferromagnetically ordered (6 = 0) ribbon in 
vacuum has a QH mean-field bandstructure [34, as shown 
in Fig. [ 1 ) 3 . (Details on the modelling are given in Ap- 
pendixjAJ. The ZLL is spin-split into the two spin sectors 
in the direction of the ferromagnetic order, denoted by 
|t), |i) at energies ±Azll/ 2 respect to the Dirac point. 
For energies within this gap, a single pair of gapless coun- 
terpropagating spin-polarised edge states develop, shown 
in red for states at the upper vacuum edge of the sample. 
Note the peculiar situation created in this energy win¬ 
dow: edge states are not chiral like in the conventional 
QH regime, but may rather propagate in both directions 
with opposite spins, like in the QSH regime. Also, val¬ 
ley degeneracy is lifted at any given edge, which hosts a 
single state per propagating direction. Upon contacting 
one edge to a conventional superconductor of gap Asc 
(as in the sketch of Fig. EF while keeping the chemi¬ 
cal potential (dotted line) within the ZLL gap, the edge 
states along the interface develop an induced gap Ag C 
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(see the corresponding Nambu bandstructure of Fig. [TJ^ 
- red lines, once more, indicate states localised at the up¬ 
per edge of the graphene ribbon, including now the dense 
quasiparticle spectrum of the superconductor). The gap 
Ag C is an important scale in this problem, since it turns 
out to be a topologically non-trivial gap. This is con¬ 
firmed by computing the bandstructure’s Z 2 topological 
invariant, Eq. (B2), relevant for quasi-one dimensional 
D -class systems [35], Unlike in a conventional QSH sys¬ 
tem, where time-reversal symmetry is required, the vac¬ 
uum edge states do not immediately develop a topological 
gap when contacted to the supercoductor, but requires 
a reasonably good contact instead. On the other hand, 
while the conventional QSH metal becomes destroyed by 
any time-reversal-breaking perturbation (such as inelas¬ 
tic scattering or magnetic impurities), which in turn spoil 
the non-trivial superconducting gap, this is not the case 
of the the present implementation, which has a broken 
time-reversal symmetry from the start. Moreover, we 
emphasize once more that no spin-orbit coupling at all is 
necessary for Ag C to develop. 

The immediate consequence of a non-trivial gap topol¬ 
ogy is the appearance of zero energy Majorana bound 
states (MBSs) at an interface with a trivial insulator, fol¬ 
lowing the bulk-boundary correspondence principle. In 
this case, however, both ends of the topologically gapped 
superconducting interface are coupled to a gapless vac¬ 
uum edge, so that the zero modes become delocalised 
into the continuum away from the interface (see the zero 
energy local density of states [LDOS] in blue in Fig. [lji). 
This situation is similar to the fate of MBSs at the ends of 
a topological proximised semiconductor nanowires when 
strongly coupled to a metallic environment. j36] 


The electronic structure associated to an antiferromag¬ 
netic ribbon (6 = 7 r, Fig. Eb is the opposite. The 
states along a vacuum edge are now (trivially) gapped 
m like the ZLL itself (Fig. [l^). Surprisingly, when con¬ 
tacting the edge to a conventional superconductor two 
pairs of gapless helical edge modes emerge with spin- 
momentum locking around conjugate momenta K and 
if* (Fig. [TJ). These unexpected states, spatially spread 
along the interface (see the blue LDOS in Fig. 01 ), are 
decoupled electron-hole (e-h) superpositions with orthog¬ 
onal and well defined spin orientation along the AF axis, 
4 ) and = a'|<^> + V\</> ht ). 

A full discussion of these states is presented in Appendix 
|Dl The two helical edge modes remain gapless as long 
as no AF canting is present in graphene (0 = n) and 
intervalley scattering is zero at the interface. We next 
consider deviations from these two assumptions. 

Canting of the AF order may be induced by means 
of a large enough in-plane Zeeman field, and is thus to 
some extent externally tunable. This idea was employed 
in Ref. [l9]to tune a graphene QH bar in vacuum between 
the AF and F regimes, leading to a insulator-to-helical 
metal transition in edge transport (evolution from Figs, 
[l^ to [TJd) . A typical canted AF bandstructure (0 = 7t/2) 
is shown in Fig. [TJi. The vacuum edge states exhibit a 


topologically trivial and ^-dependent gap, smaller than 
the bulk Azll- Along a superconductor interface, the 
canted AF helical states are also gapped, Fig. [l}. Like 
in the ferromagnetic case, this gap is topologically non¬ 
trivial. This situation allows for the emergence of true 
localised zero-energy MBSs at the ends of the supercon¬ 
ductor interface, where the edge gap changes topology, 
see Fig. IT The MBSs are topologically protected, and 
are not destroyed by any small perturbations, or even 
by modifying the crystal structure of the superconductor 
(Appendix |E|). 

To understand the full phase diagram of the 
graphene/superconductor interface quantitatively, it is 
useful to employ a simplified description in terms of an 
effective low-energy Hamiltonian for the edge states (see 
Appendix |B|. The model is valid for a chemical poten¬ 
tial tuned to the ZLL gap, and has the advantage of 
allowing us to incorporate the effects of atomic disor¬ 
der along the junction (encoded in an intervalley cou¬ 
pling w, where ‘valley’ here refers to the conjugate if 
and if* momenta) and arbitrary spin canting (encoded 
in an intravalley splitting be = Ag C cos[0/2]). It correctly 
describes the three possible phases for low-energy inter¬ 
face modes: gapless, trivially gapped and non-trivially 
gapped. The corresponding phase diagram is shown in 
Fig. [Ip. Typical edge-mode dispersions within each 
phase (with if and if* points folded onto the T point) 
are shown in Fig. [TJ , and are characterised by their en¬ 
ergies /a i ? 2 < Ag C at k x = 0 for 0 = tt [AF], w = 0, and 
their corresponding velocities v\^ >0 (left panel). 

The gapless interface regime (light blue in Fig. If) is 
achieved for \w\ < Wi ns = ^ l/iiyAq/tq — l^ 2 \/W/v 2 \ and 
6 > 6> ins , where b 6ins = A£ c cos(6> ins /2) = \\fiiyjv 2 /vx + 
(12 /V 2 1 • The gapped regimes are characterised by 

the Z 2 topological invariant of the system, which reads 
v = sign (w 2 + /ii/i 2 — tig) (see Appendix |b| . A trivially 
gapped phase v = +1 is reached for strong intervalley 
coupling w at the interface, w hile for int ervalley scatter¬ 
ing below a threshold w < yjb 2 e — /ii/i 2 , the interface is 
one-dimensional TS with invariant v = — 1. The non- 
trivially gapped regime is most robust against disorder 
for F order, for which the threshold w reaches its maxi¬ 
mum wq = yj (Ag C ) 2 — /ii/i 2 . Note that to achieve a non¬ 
trivial TS interface, the intervalley coupling w should 
therefore never exceed the induced gap Ag C (this is al¬ 
ways the case for sufficiently transparent junctions). 


III. EXPERIMENTAL SIGNATURES OF 
GRAPHENE MAJORANAS 

We finally consider measurable signatures of the MBSs 
in the system. A powerful probe whose feasibility has 
been recently demonstrated experimentally [ 20 | [ 21 ] in¬ 
volves interferometry of critical currents in Josephson 
junctions, and Fraunhofer pattern anomalies in partic¬ 
ular mmm • Lee et al. predicted [38] that topological 
superconductivity in a short and wide Josephson junction 
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FIG. 2. Normalised critical current I c /Ic as a function of 
magnetic flux through a Josephson junction. Curves from 
bottom to top (shifted for better visibility) correspond to the 
non-interacting case (black), gapless AF phase (purple), a 
trivially gapped AF phase (light gray) and a ferromagnetic 
phase with a topologically non-trivial gap (orange). Only the 
latter shows non-decaying non-zero minima in / c , a conse¬ 
quence of an Andreev spectrum (inset a) with an odd num¬ 
ber (one) of edge-resolved zero energy crossings (bottom-right 
panel). States coloured in red and blue are located at the top 
and bottom edges, respectively (inset b ), while states in gray 
are spread across the width of the junction (inset c). 


could be directly detected in its Fraunhofer pattern, in 
the form of non-vanishing minima of the critical current 
for arbitrary magnetic flux through the junction. Such 
Fraunhofer anomaly was recently observed in a three- 
dimensional topological insulator, and was interpreted as 
possible evidence of MBSs [12]. 

Fig. |5| shows the Fraunhofer pattern in a short 
and wide (lOnm x 3/im) graphene Josephson junction 
in various regimes. The black (bottom) curve corre¬ 
sponds to the non-interacting case (no magnetic or¬ 
dering, Azll = 0), which exhibits the conventional 
J c (4>) = J°| sin(7r4>/4>o)|/(7r4>/4> 0 ) critical current that 
decays as the inverse magnetic flux through the junc¬ 
tion I c ~ l/4> and vanishes at multiples of the flux quan¬ 
tum 4 >q = h/2e. A similar behaviour is obtained in the 
gapless regime 0 > 0[ nS: w < w- ins (purple curve, with 
0 = 7r [AF] and w = 0). In both cases, the junction 
is host to a narrow quasi-continuum of Andreev bound 


states around the Fermi energy for any value of the super¬ 
conducting phase difference <fi. The corresponding spec¬ 
tra are shown in the top row of inset a (for = 15.54>o). 
States in red and blue are localised at the top and bot¬ 
tom edges of the junction (inset 6), respectively, with 
gray denoting states spread across the junction (inset c). 
The case with trivially gapped interfaces (strong interval¬ 
ley scattering) also exhibits a generic Fraunhofer pattern 
with vanishing minima (light gray curve, 6 = 0 [F] and 
w > wo)- The minima, however, occur at fluxes that 
are shifted away from integer 4>/4>o at high <F, while the 
maxima do not decay like the conventional I c ~ l/4> pat¬ 
tern, which is connected to non-uniform currents across 
the junction 37 . The Andreev spectrum of the trivially 
gapped phase is qualitatively different from the gapless 
spectra, and generally shows a distinct gap devoid of any 
edge states (inset a, bottom left). For certain values of 
parameters, it may exhibit zero-energy crossings inside 
the gap, but in such cases these crossings are accidental 
(not topologically protected) and there is always an even 
number of them at a given edge. 

The Josephson junction with a non-trivial gap along 
the contacts is distinctly different from all previous cases. 
This phase develops two MBSs at each edge (top and 
bottom) which hybridise to carry a finite supercurrent 
that never vanishes as the flux increases. The corre¬ 
sponding Fraunhofer pattern thus exhibits a finite back¬ 
ground with a superimposed non-decaying oscillation (or¬ 
ange curve in Fig. [2]), as described in Ref. [381 The 
finite minima are roughly one half of the maxima at 
large flux, and occur away from integer 4>/4>o, at val¬ 
ues very close to the zeroes of the trivially gapped phase. 
(Note, however, that such distinctive pattern develops 
only in junctions shorter than the Majorana localization 
length and pierced by a large number of flux quanta, 
so that the Major anas are well developed and opposite 
edges are decoupled). These Fraunhofer anomalies, al¬ 
though not completely unambiguous, thus constitute a 
measurable hint of the presence of MBSs at the end of a 
graphene/superconductor interface. Unfortunately, fab¬ 
ricating a very short junction is challenging, in particular 
due to charge-transfer effects from the superconductors 
which were neglected here, and which will dope graphene 
away from neutrality within a few nanometers of the con¬ 
tacts. It is thus important to explore other less stringent 
experimental schemes that are at the same time not am¬ 
biguous. The key is to probe the Andreev spectrum di¬ 
rectly for signatures of Majoranas and non-trivial topol¬ 
ogy* 

The presence of the two hybridised Majoranas per vac¬ 
uum edge in the non-trivial phase manifests in the An¬ 
dreev spectrum as a single topologically protected zero 
energy crossing at each edge as </> is increased by 2tt (one 
red and one blue crossing, see bottom-right inset a in 
Fig. §. An odd number of such zero energy crossings 
has been shown |39j to be an direct manifestation of non¬ 
trivial topological order v =» — 1, and is the underlying 
reason for the anomalous Fraunhofer pattern of the junc- 
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FIG. 3. (a) In blue, density of the lowest Andreev level formed by the hybridization of four Major ana bound states in a 
square 500 nm x 500 nm Josephson junction at different values of canting angle 0. Rest of parameters like in Fig. ^ In A 
and C (0 < 0l and 0 > 0w respectively) the hybridization is strong, while in B it is exponentially suppressed, and the level 
remains a four-fold degenerate zero mode, (b-e) Transport spectroscopy dl/dV from a normal point contact N (transparent, 
single spinful channel), see (a), (b) and (d) show (at 0 — 0 and 0 = 0.17T, respectively) the dl/dV as a function of junction 
phase difference 0. Both exhibit an odd number (one) of zero energy crossings of edge-resolved dl/dV resonances, signalling a 
topologically non-trivial system, (c) dl/dV as a function of canting angle 0 at 0 = 0. It exhibits a zero bias anomaly in the 
window 6l < 0 < 0w , see dashed lines, (e) The dl/dV dependence with out of plane magnetic field B z through the sample 
(0 = 0, 0 = 0.l7r) shows the same zero-crossings as with 0 each time the total flux <f> = B Z LW is increased by a flux quantum 
Mo¬ 


tion. A completely non-ambiguous demonstration of the 
presence of MBSs is also thus possible in principle, by di¬ 
rectly counting edge-resolved zero-energy crossings using 
Andreev spectroscopy m in a phase-controlled Joseph¬ 
son junction. This may be achieved by measuring differ¬ 
ential conductance dl/dV through a normal point con¬ 
tact attached to one edge of the junction, as sketched in 
Fig. I a). We assume a single spinful channel is open. 
Each Andreev level of energy e in the junction is detected 
as a dl/dV resonance through the probe at bias V = e/e, 
with a resonance width that measures the state’s proba¬ 
bility density at the point contact. Figs. ,d) show a 
simulation (see Appendix [C] for details) of such a dl/dV 
as a function of 0 (at 0 = 0 and 6 = O.ln respectively) 
for a square 500 nmx 500 nm Josephson junction. Note 
that, unlike in the Fraunhofer simulation, this is not a 
short junction, since that is no longer a desirable or real¬ 
istic requirement in the context of Andreev spectroscopy. 
The number of edge-resolved zero energy crossings as 0 
is swept from 0 to 2tt is one, as corresponds to Majorana- 
hosting SC contacts (compare this dI/dV{4>) to the blue 
lines in bottom-right inset Fig. |2^i [short junction]). Inci¬ 
dentally, the dl/dV (0) profiles at zero (non-zero) 0, panel 
b (d), follow the characteristic Andreev level spectra in 
topological Josephson junctions through long semicon¬ 
ducting nanowires at perfect (non-perfect) transparency 
[4TH45]. The conductance at the crossings is pinned to a 


universal value 4 e 2 /h (white spot), see Appendix [cj 

While a 0-controlled junction typically requires a 
SQUID-like geometry and may be experimentally chal¬ 
lenging, the existence of MBSs in the junction may be 
detected even more simply by varying the out of plane 
magnetic field B z , and hence the total flux through 
graphene. An increase of by 4>o is equivalent to in¬ 
creasing 0 by 27r for large This is shown in Fig. 

§ (compare to Fig. § 0 - Moreover, at fixed B z and 
0 = 0, the MBSs also show up as a zero-bias dl/dV peak 
(Fig. [ 3 J 3 ) as 0 is tuned - by the in-plane Zeeman field 
B\\ - within a range 0l < 0 < #w < #ins (dashed lines). 
For the realistic parameters used in Fig. [ 3 ] (B z = IT, 
L = W = 500nm), and using Bq = 10T as the typ¬ 
ical in-plane field B\\ for complete Ferro polarization, 
6l ~ 0.177T and 0w ~ 0.57T are reached for Bu ~ 9.6T 
and B\\ « 7.IT, respectively (see Appendix 0. Inside 
this window (point B in Figs. pm,c), the four Majoranas 
are concentrated at the corners of the junction and do 
not overlap, as they decay within a distance smaller than 
both the width W and length L of the junction, and hence 
appear as a sharp zero-bias resonance in the dl/dV with 
exponentially small splitting. The resonance has a uni¬ 
versal magnitude of 2 e 2 /h at low temperatures, see Ap¬ 
pendix^ This is analogous to the zero-bias anomalies re¬ 
ported in pioneering Majorana experiments on semicon¬ 
ducting nanowires pQ. In contrast, for 0 < Ol (point A, 
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bulk approaching ferro-ordering), Majorana pairs over¬ 
lap along the vacuum edge and produce a split resonance. 
Likewise for 0 > 6w (e.g. point C, SC contact close or in¬ 
side the gapless regime), Majoranas overlap along the SC 
contact, and develop a (roughly (^-independent) splitting, 
making the 0 > Ol junction strictly trivial. Note also the 
strong suppression of the width in the dl/dV resonances 
in this case, due to the exponentially small wavefunction 
amplitude at the point contact in this geometry (dotted 
lines overlaid for 6 > 0.67T to improve visibility). Essen¬ 
tially the same dl / dV ( 0 ) phenomenology is obtained in a 
setup with a single superconducting contact (which hosts 
two Majoranas instead of four), see Appendix [Cj 


IV. DISCUSSION 

Our results show that the spontaneous magnetic or¬ 
dering of the ZLL in graphene enables the creation of 
topological superconductivity and Majorana states at an 
interface with a conventional superconductor, even in the 
absence of spin-orbit coupling in the system. The key is 
to tune the Fermi energy in the contact into the ZLL 
gap, and to achieve a good proximity effect therein. The 
recently characterised samples of Ref. [20] are good can¬ 
didates to realise our proposal. Impressive progress in 
controlling graphene filling into proximity gaps has also 
been reported [46] . We furthermore showed that non¬ 
vanishing and non-decaying supercurrent minima in the 
Fraunhofer pattern across a depleted graphene Joseph- 
son junction constitute a characteristic signal of topolog¬ 
ical order and the presence of Majorana bound states in 
the junction [38]. Fraunhofer patterns of extraordinary 
quality have been recently reported in high-transparency 
ballistic graphene Josephson junctions [21]. We predict 
even stronger observable signatures of non-trivial topol¬ 
ogy in Andreev transport spectroscopy, both in the form 
of an odd number of 4 e 2 /h edge-resolved zero-bias cross¬ 
ings versus junction phase difference </> or out of plane 
magnetic field L>, and extended 2e 2 /h zero bias anoma¬ 
lies versus canting angle 0. These experimental probes 
and the required device parameters are within reach in 
top laboratories today. We thus expect that the possibil¬ 
ity of tuning graphene/superconducting interfaces into a 
topological phase hosting Majorana bound states could 
be tested soon. 
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Appendix A: Modelling 


In this section we present the system models employed 
in this work. A non-interacting graphene flake may be 
modelled by a nearest-neighbour tight-binding Hamilto¬ 
nian in an honeycomb lattice, with lattice constant ao 

#0 = - F tel4>ijc ]s c is + N m n is (Al) 

(i,j),s is 

where rii s = cJ s Ci 3 , i is the site index, s is spin, and 
ct>ij = — | dr • A(r) is the Peierls phase due to the 
magnetic flux B z = z • (V x A). We consider a perturba¬ 
tion JT B • S] arising from an external Zeeman field H, 
where Si = ^2 ss f c\ s a SS 'Ci S / is the spin at site i. Intrin¬ 
sic electron-electron interactions is furthermore included 
in the local Hartree-Fock approximation eh. which then 
take the form of a self-consistent Zeeman-like field Bjj{ri) 
that is different in the two honeycomb sublattices. The 
total Zeeman-like perturbation Hz thus reads 


H z = Yl [B + BuFii 

i 


Si 


(A2) 


While B is uniform, favoring ferromagnetic (F) ordering, 
the self-consistent Bu{fi) is generally opposite for nearest 
neighbours, favouring antiferromagnetic ordering of the 
bulk. The combination of the two leads to a H-tuneable, 
spin-ordered ZLL that can be tuned from AF to F, as 
discussed in the main text. Further details on the mean 
field numerics and results can be found in Ref. EU 
The hybrid graphene/superconductor (SC) system is 
described by the Hamiltonian 


H = H 0 + H z + ^2 [AsccjfT + h - c 

r=i€SC 


(A3) 


The Fermi energy in H 0 above is (i^ for r* ^ SC in 
graphene and (is for fi G SC at the superconductor (also 
a honeycomb lattice here, although this is not essential, 
see Appendix [e|. Similarly L>, Bu and B z are zero in 
the superconductor. However, gauge invariance demands 
that for a finite magnetic flux in graphene, A seif) = 

A sc ex P — f r df’A(r ) , where A sc is the pairing for 


zero flux at a given superconductor. 

The critical currents for the Fraunhofer patterns have 
been calculated in a wide and short graphene Josephson 
junction, described by a discretized H using an up-scaled 
ao for numerical efficiency, following the ideas of Ref. [371 
The critical current for each magnetic flux is calculated 
as I c = 2e/fix ma x^dF/d(/)), where (j> is the supercon¬ 
ducting phase difference and F(</>) is the free energy of 
the junction. [48] Exact diagonalization of the Hamil¬ 
tonian is used to evaluate F(0) at zero temperature. 
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The method to compute the differential conductance for 
transport spectroscopy is explained in Appendix |C] 


Appendix B: Low-energy description of a 
graphene/superconductor junction 


In this Appendix we derive a simplified effective Hamil¬ 
tonian for the two helical edge modes below the su¬ 
perconducting gap Asc that arise in a generic Quantum 
Hall graphene ribbon and a superconductor. We also 
characterise its topology by deriving expressions for the 
relevant topological invariant as a function of model pa¬ 
rameters. 

We assume the chemical potential lies within the gap 
Azll induced by interactions in the zero Landau Level 
(ZLL). The gap is associated to a bulk spin ordering 
described by a canting angle 6 between the two sublat¬ 
tice. The effective model incorporates an arbitrary value 
for 0 in graphene and also intervalley coupling due to 
atomic disorder along the interface. Formally, H e ^ is a 
projection of the microscopic Hamiltonian on the basis 
{\Kt), | Ki), \K*l)} of the four AF helical edge 

states (see Fig. If in the main text, and their analytical 
description in the preceding section) along the junction. 
We furthermore consider a linearisation of their disper¬ 
sion in the AF case around ‘valleys’ K and K* points. 
These two valleys are folded onto the Y point by appropri¬ 
ately expanding the ribbon unit cell. Such folding allows 
us to include intervalley scattering into i7 e ff in a simple 
way. H e ff then takes the form (ft as 1) 


1. Topological invariant of edge Hamiltonian 


In symmetry class D (superconductors without time 
reversal symmetry), the one-dimensional topological in¬ 
variant is Z 2 . It is conventionally defined, for periodic 
systems, as [49] 


. (Yf[H(0)r x ] 


So_ 

S^TT 


(B2) 


where s a = signPf [H(a)r x \, H(k x ao ) is the ID Bloch 
Hamiltonian for momentum k x , ao is the lattice con¬ 
stant of the ribbon lattice, r x is the first Pauli matrix 
in the electron-hole sector, and Pf is the Pfaffian. It 
can be shown in general that Hr x is antisymmetric at 
the high-symmetry points k x ao = 0,7r. The invariant v 
is thus fully determined by the structure of H at these 
two points. The effective Hamiltonian i7 e ff 5 Eq. (Bl), 


only gives a faithful representation of the full microscopic 
Hamiltonian H around one of them, the folded Y point 
k x = 0 (Fig. lj, main text). To extract analytic results 
for v for the full H using i7 e ff, we must ensure that at 
the 7r-point s n does not change when sweeping the pa¬ 
rameter space (since this sector of states is not described 
by FZeff). This is indeed the case in our system for the 
chosen basis, for which s n = 1. The changes in topology 
stem from the reconnections of the low-energy edge states 
that are concentrated around T, and are well described 
by i7 e ff- Higher excited states not included in H e q never 
cross zero energy, and therefore cannot affect the sign of 
the Pfaffian of H(0)r x . One can thus write 


v = signPf[iJ e ff(0)r a: ] = sign(V + \i hj, 2 - b 2 ) (B3) 


M 1+^1 k x bg 


w 


M 2 - V 2 k x 0 w 

0 -M2 - v 2 k x bg 

w bg -M 1 +V 1 k x 

(Bl) 

ey coupling bg = Ag c cos(6>/2) 
implements AF canting, and couples opposite spins 
within the same valley. The intervalley coupling w is 
spin-independent, and corresponds to the harmonic of 
wavenumber A K = ( K * — K) • x of any disorder term W 
close to the interface, w = ((f>^K*\W(AK)\(/)^K}- Both 
bg and w can be chosen real without loss of generality. 
^ 1,2 > 0 are the velocities of the counter-propagating he¬ 
lical states, and /ip 2 < Ag C are their energy, relative 
to the Fermi energy, at the Y point (see Fig. lj in the 
main text). The overall structure apparent in H e q is fully 
determined by the particle-hole symmetry of the under¬ 
lying Nambu description. Note that H e ^ only retains 
terms linear in momentum k x along the interface. 

The simplicity of the linearised low energy model above 
allows us to compute analytical expressions for the topo¬ 
logical invariant and the boundaries that separate dif¬ 
ferent phases of the junction (gapless, trivially gapped, 
non-trivially gapped), as summarised in the main text. 
We now briefly sketch their derivation. 


tfeff 


be 

w 

0 


Here, the intravail 


We have numerically verified the above result by evaluat¬ 
ing the Z 2 invariant exactly from the microscopic Hamil¬ 
tonian H. 

We finally sketch the derivation of the insulating 
thresholds u;i ns and 0- ins . These are extracted by comput¬ 
ing the solutions for the wavenumber k x of H e ff modes at 
zero energy. Since H e g(k x ) is linear in k x , said k x solu¬ 
tions at e = 0 can be obtained as eigenvalues of a matrix 
— (9/c x i7eff(0)) _1 i7eff (0). These can be worked out ana¬ 
lytically, and turn out to be all complex (i.e. the interface 
becomes gapped, see e.g. Fig. [7]f) if w > W[ ns or 6 < #i ns , 
with the expressions given in the main text, 


^ins — ~ 


Mi - M2 y/v-i/v 2 


(B4) 


= A sc cos(6» ins /2) = t Mi y/v 2 /vi + M 2 ^ 1/^2 


Appendix C: Transport spectroscopy 

1. Computing differential conductance 

The computation of the transport spectroscopy results 
shown in Fig. [3] follows standard techniques of quan¬ 
tum transport. All interactions are incorporated into the 














FIG. 4. (a) 350 nm x 350 nm sample with a single superconducting contact, otherwise identical to system of Fig. [3] main 
text, (b) and (c) show the differential conductance as a function of canting angle 0 and out of plane magnetic field B z . Both 
show the same phenomenology as panels (b) and (e) of the two contact setup in Fig. [3j albeit with a simpler two-Majorana 
geometry, and with the total length of the vacuum edge L* = 2L + W playing the role of L. 


mean field solution for Bu , Eq. ( |A2| . In this case, the 
differential conductance dl/dV through a normal con¬ 
tact with a single spinful channel at a certain bias V is 
given by the BTK formula dl/dV = (2 — R ee -\- Rh e )e 2 /h, 
where the total electron-electron (R ee ) and electron-hole 
(Rhe) reflection probabilities of free electrons incident on 
the contact are evaluated at energy e = eV. The full R 
matrix, including both electron and hole sectors, may be 
obtained in terms of Caroli’s formula R = Tr (G r rG a r), 
where G r / a are the (dressed) retarded/advanced Green 
functions in the sample, and T = — i(E — E*) is (twice) the 
decay rate matrix into the normal probe. All these op¬ 
erators are defined in the Nambu basis, just like e.g. Eq. 
(D5). The Green’s function matrices where solved by in¬ 
verting the Dyson equation (e ±i0 + — H — E)G r / a (e) = I 
using efficient linear algebra routines, where H is the 
Nambu Hamiltonian of graphene, including the supercon¬ 
ductors. The self-energy E = V^g r V, defined in terms of 
the hopping matrix V between graphene and the semi¬ 
infinite normal lead, and the retarded Green function g r 
of the latter. g r is obtained by solving the corresponding 
(self-consistent) Dyson equation (e + i0 + — h — v^g r v)g r = 
I, where h and v are the on-site and hopping matrices 
acting on the lead’s constituent unit cells. 


Fig. IJa. Likewise, 0w is defined as the 6 such that the 
corresponding decay length of edge states along a super¬ 
conducting edge equals the width W (hence 0w < #ins)- 
The evaluation of 0w can be made by extracting the 
decay length l/lmk x of gap state s at zero energy from 
the effective Hamiltonian of Eq. (Bl) without disorder 
{w = 0), and equating that to W. The result comes out 
simply as 


cos 0w ~ cos #i ns + 


2viv 2 
W 2 A* 2 C 


where, recall, h= 1 and 


COs(#j ns /2) = i fj,l^/v 2 /v 1 + Il 2 \/vijv 2 /A; 


sc 


Note that in the limit W oo, 0w = ^ins, as expected. 

An analogous calculation can be done for the vac¬ 
uum edge along the y direction, whose effective (no rmal) 

h ® 


Hamiltonian can be written in analogy to Eq. 


as 


irvac 
-Cf eff 




flN + Vpky 

zll sin(#/2) 


^Azll sin(6 > /2) 
Mat — vpky 


(Cl) 


This leads to the estimate 


2. Estimates for canting angles 6 


L,W 


■ 2 0l _ / 0 Mat 
sm — ~ ( 2 


*ZLL 


2vp 

LA 


ZLL 


The canting angle 0l is defined as the 6 such that the 
corresponding decay length of edge states along a vacuum 
edge equals the length L of the Josephson junction, see 


Note that as L oo, 0l reaches a minimum value that 
corresponds to the threshold where the vacuum edge be¬ 
comes gapped (non-zero for ^ 0). 
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FIG. 5. Transport spectroscopy dl/dV obtained in various arrangements of normal, trivial superconducting (SC), and topo¬ 
logical superconducting (TS) nanowires, to be compared to Fig. [3j The wires are modelled after Refs. MandflTl (a) A trivial 
SC/normal/SC Josephson junction, is probed by a third normal contact at bias V, as a function of superconducting phase 
difference (/). No protected zero-bias dl/dV anomalies aries. (b) Same as (a) for a TS/normal/TS junction. A universal 4 e 2 /h 
zero bias anomaly (white) is obtained at a certain 0, here (j) — 7r, for which the two Majoranas in the junction (blue circles 
in the inset) decouple, (c) A normal/SC junction that transitions into normal/TS with two Majoranas as the longitudinal 
Zeeman field Vz exceeds a critical value V%- Apart from a small oscillatory splitting that decays exponentially with TS length, 
the dl/dV in the latter case shows a universal 2 e 2 /h zero-bias anomaly (orange). 


Mean field results within the Hubbard model (see Ap¬ 
pendix A] and Ref. ED yield a dependence of cant¬ 
ing angle 6 with in-plane magnetic field B\\ of the form 
sin(0/2) « yjl — (T>||/T>o) 2 , where B 0 is the in-plane field 
that achieves complete Ferromagnetic polarization. 


3. Differential conductance with a single 
superconducting contact 

The formation and detection of graphene-based Majo¬ 
ranas only requires a single superconducting contact. In 
that sense, the geometry discussed in Fig. [3] of the main 
text, while relevant in the context of the Josephson effect, 
Fraunhofer patterns and parity crossings, is not minimal. 
A simpler geometry with a single superconducting con¬ 
tact allows for the detection of two Majoranas (instead 
of four) in the form of a zero-bias anomaly, analogous 
to that of Fig. [3^. In Fig. [4] we present the dl/dV 
in such a geometry. Panel b shows the formation of a 
zero bias 2e 2 /ff anomaly within 0\y < 0 < Ol* (with a 
L* = 2L + W that now corresponds to the total length 
of the vacuum edge). Note that, while this setup (panel 
a) does not allow for a (^-controlled modulation, its 0 
dependence is qualitative the same as for a Josephson 
junction. Interestingly, moreover, the analogous to the </> 
modulation induced by changing the flux B z does oper¬ 
ate in this setup just like in a Josephson junction. The 
reason is that the flux changes the relative phase of the 
two Majoranas in the SC contact, making them cross at 


zero energy each time the flux is increased by ff>o. This 
is shown in Fig. II- 


4. Relation of the graphene dl/dV to junctions of 
topological nanowires 

The transport spectroscopy results for Majoranas in 
graphene presented in the main text exhibit three differ¬ 
ent regimes, labeled as A, B and C in Figs. [3] and [4] 
These arise due to the interplay, as a function of cant¬ 
ing angle 0 between delocalization of Majorana bound 
states along either a vacuum edge and the supercon¬ 
ducting contact. From the point of view of the normal 
point contact, the former case (A) is analogous to a one¬ 
dimensional TS/normal/TS Josephson junction, where 
the normal probe is an extra lead coupled to the normal 
section for spectrocopy. In contrast, in the case for which 
the Majoranas remain bound to the corners of the sam¬ 
ple and do not delocalize (R), the probe is tunnel coupled 
to the closest Majorana bound state (top-left corner of 
the sample), and plays the role of a tunnel normal/TS 
junction like in the zero-bias anomaly experiment of Ref. 
[T|. Finally, case C is like case B albeit for a trivial nor- 
mal/S junction without Majorana bound states. These 
mappings to well-understood systems are useful to un¬ 
derstand the universal values of the dl/dV obtained in 
each case. 

For case analogous toff, a one-dimensional normal/TS 
junction, it is well known [50U56] that the dl/dV for a 
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long enough TS yields a universal 2 e 2 /h zero-bias con¬ 
ductance resonance, a telltale signal of the presence of a 
Majorana bound state at the contact. If the TS has a 
finite length, the overlap of the Majorana with its sibling 
at the opposite end of the TS section gives rise to a split¬ 
ting of the zero-bias resonance by an energy exponentially 
small in the TS length divided by the spin-orbit length 
(or, more precisely, the Majorana localization lengthj57]). 
A simple model for such a N/TS junctions was devised 
by Oreg et al. and by Lutchyn et al. in Refs. 16' and □n 
The model is based on semiconducting wires that exhibit 
a TS phase when proximized to an s-wave SC while un¬ 
der a longitudinal Zeeman Vz exceeding a critical value 
Vg = yj/i N + A 2 . The typical dl/dV in such a junction 
as a function of Vz and bias V is shown in Fig. Hf- Note 
that, indeed, for Vz > V%, a zero bias anomaly of mag¬ 
nitude 2 e 2 /h develops, with a small splitting (oscillating 
in Vz) EBM]. 

In the case analogous to A, we have a TS/normal/TS, 
where the normal portion represents the graphene vac¬ 
uum edge over which the Majoranas delocalize, and 
probed with an additional point contact. The Oreg- 
Lutchyn model in the topological phase yield a trans¬ 
port spectroscopy map, shown in Fig. with a single 
zero energy crossing as the junction phase difference </> 
increases by 27r, just like in Fig. !b ,d), by virtue of 
the non-trivial junction topology. Moreover, the zero- 
bias dl/dV at the crossing is pinned to 4 e 2 /h (white), 
again like in the graphene case Fig. [3](b,d) . This univer¬ 
sal value can be understood intuitively as the addition 
in parallel of two 2 e 2 /h normal/TS zero bias anomalies, 
one per Majorana (both are coupled to the probe in this 
geometry), when the two become decoupled at the ap¬ 
propriate (j). 

Finally, note that the trivial S/normal/S Josephson 
junction, corresponding to case (7, does not yield a uni¬ 
versal zero-bias anomaly at any phase 0, see Fig. [5^t. In¬ 
stead, the finite energy Andreev levels yield a dl/dV at 
finite bias that approaches (non-universal) 4 e 2 /h. Un¬ 
like for the topological case A above, however, pertur¬ 
bations to the system may introduce additional normal- 
reflection component to said Andreev levels that suppress 
this value. 


Appendix D: Helical edge states in an AF 
graphene/superconducting interface 


1. Interface states without Landau levels 

The interface states between a superconductor and an 
antiferromagnetic honeycomb lattice are not related to 
the Landau level structure. In the particular case of 
graphene, the magnetic field is the key ingredient to de¬ 
velop magnetic order (due to the large kinetic energy of 
electrons), which is developed when the kinetic energy is 
quenched by the magnetic field. However, in a general 
honeycomb lattice, provided the interactions are large 


enough at B = 0, electrons might be able to develop AF 
order, and thus create interface states between a SC even 
at B = 0. This behaviour might be relevant for other 
honeycomb systems with smaller hopping strengths than 
graphene, as silicene, germanene, stanene or honeycomb 
oxides. 

We show in Fig. [6j that such interface sustains the 
same kind of states as in the case of the antiferromag¬ 
netic quantum Hall state. In the case of zigzag inter¬ 
faces ©#)> each valley supports its own set of inter¬ 
face states, whereas for an armchair interface ©>), the 
two valley are folded and interface states between differ¬ 
ent valleys can couple. In the former case, if the interface 
is abrupt enough, a small gap opens up due to intervalley 
mixing. 

If a canting in the magnetic moments is introduced, 
the zizgag interphase remains gapless. The same happens 
when only a orbital magnetic field is introduced. Only 
when both perturbations are present simultaneously, the 
system is able to enter into the topological superconduct¬ 
ing state. Thus, an off-plane magnetic field is mandatory 
to observe the Majorana bound states. The previous phe¬ 
nomenology, suggests that in order to develop a topologi¬ 
cal gap, both the spin rotation symmetry and the spatial 
gauge symmetries have to be broken. 


2. Helical edge states from wavematching 

The interface states between a honeycomb antiferro- 
magnet and a superconductor are not intrinsically related 
to the Landau level spectrum. Although in graphene, the 
antiferromagnetic state is only expected to arise when 
the system enters in the quantum Hall regime, a general 
antiferromagnetic honeycomb lattice might also sustain 
interface states when attached to a superconductor with¬ 
out a magnetic flux. 

In this section, we will show how that interface states 
naturally arise by an analytic argument in the absence 
of magnetic field. In particular, a simple wavematching 
between a E = 0 energy state shows that the boundary 
between an antiferromagnet and an swave superconduc¬ 
tor is able to sustain such state. 

To proceed, we will look for bounded solutions such 
that H\4>) =0 with 


(t>(x) = 


r(x) = 


0 -\(x)x 


(Dl) 


X(x) = 



if x < 0 
if x > 0 


(D2) 


In the following we will focus in one of the four de¬ 
coupled sectors, in particular the |e, t K) with |ft, 4- K') 
sector 
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FIG. 6. (a) Scheme of a ribbon hybrid ribbon AF-SC, with 

the color code of the bandstructure. (b) Band structure of a 
hybrid ribbon with doped SC and armchair interface. Band 
structure of a hybrid ribbon with un-doped (c) and doped (d) 
SC with zigzag interface, which shows the different interface 
states in each valley. 



0.4 

0.2 

w 0 

- 0.2 

-0.4 


4>(x) 


1 \ 

i 


(D7) 


e Ascx if x < 0 

e~ rnx if x > 0 


(D 8 ) 


so that a normalizable E = 0 exist for an interface 
between a trivial antiferromagnet and a trivial Dirac su¬ 
perconductor. 


3. Helical modes from explicit integration 

In previous section, we build the E = 0 by wavematch¬ 
ing across a sharp interface. However, it is is possible to 
give a general solution for the interface state between 
the antiferromagnet and the superconductor. Without 
loss of generality, in the following we will assume m > 0 
and A sc > 0 . The Hamiltonian for an arbitrary antifer¬ 
romagnet and pairing profile for k y = 0 reads 


For the antiferromagnet the Hamiltonian reads 


H S c = 


/m p 0 0 

p —m 0 0 

0 0 m —p 
\ 0 0 —p — rnj 


whose E = 0 solutions are 


(D3) 


H = 7 i m(x) + 72 p + 73 A S c(x) 
with 71 , 72,73 defined by 


7i = 


/I 0 0 0 \ 
0-100 
0 0 10 
\0 0 0 - 1 / 


(D9) 


(D10) 


V2 


i 

0 

V 




(D4) 


On the other hand, for the superconductor the Hamil¬ 
tonian reads 


H S c = 


with E = 0 solutions 



(D5) 


4,1 V2 


* 2 = 7I 


(D 6 ) 


Imposing continuity at the interface x = 0, the full 
E = 0 solution reads 


72 


/0 1 0 0 
10 0 0 
0 0 0 -1 
Vo 0 -1 0 


73 = 


/0 0 1 0 \ 
0 0 0 1 
10 0 0 
Vo 1 0 0/ 


(Dll) 


(D12) 


Defining 74 = —^7271 and 75 = ^ 7273 , the zero energy 
equation reads 


(p + - i^Ascix))^ = 0 

The spinor wavefunction 


(D13) 




1 


(D14) 


V i 
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FIG. 7. Scheme of the effect of chemical doping in the Dirac 
spectrum of proximized graphene (a). Schematic phase di¬ 
agram (b) of the interface electronic spectrum between the 
SC and the QH ferromagnet, as function of the in-plane field 
and the chemical doping of the SC. Quasiparticle energies 
for a gapless (c) and gapped (d) interface. Phase boundary 
obtained by numerical calculation at low dopings (e) and at 
large doping (f). The green and red colored states of the band 
structures correspond to the states localized at the interface, 
where the topological superconducting gap is calculated be¬ 
tween the red states shown in (c,d). The blue gapless states 
correspond to the chiral states between the QH and vacuum. 


verifies 74^0 = 7500 = 0ce which allows to build the 
E — 0 solution 


a 0.4 
0.2 


- 0.2 
-0.4 

0 0.2 0.4 0.6 0.8 1 

k/(2n) 






FIG. 8. Bandstructure of the an interface between a square 
superconductor and an antiferromagnet (a) or canted antifer- 
romagnet (c) honeycomb lattice. Figures (b) and (d) show 
the local density of states in a finite system, which corre¬ 
spond to gapless interface channels (without canting) (b) and 
Majorana states (with canting) (d) 


leads to a one dimensional gapless state (see Fig. [ 7 ]^) . In 
order to reach the interfacial topological superconducting 
state, the bended bands (Fig. [7^) have to be moved up 
in energy. This can achieved by increasing the in-plane 
field, so that 0 < 0[ ns . The critical in-plane field B x 
as a function of doping, which separates the gapless and 
topologically gapped states is shown in Figs. [7^,f. For 
small Use, the critical field increases linearly, leading to 
small critical fields at small doping in the SC, whereas 
for large doping the critical field saturates. 


4>(x) = Ne~ fo( m ( x ')- A sc(x'))dx' fa ( D15 ) 

which is normalizable provided that m(+oo) > 
^S'c(+ oc ) and Asc(—°°) > m(— 00 ), which is the con¬ 
dition of domain wall between superconductor and anti¬ 
ferromagnet. For the case of step profiles, the solution 
obtained by wavematching is recovered. 

4. Influence of fisc in the critical Zeeman 

In the idealised situation in which the superconduc¬ 
tor is described as a single-orbital honeycomb lattice at 
half filling, a arbitrary small Zeeman field is capable of 
opening the interface topological gap. However, charge 
transfer processes are expected to shift the chemical po¬ 
tential of the proximized graphene (SC region in Fig. [l]). 
In this situation, band bending of the interfacial states 


5. Emergence of AF helical states from a 
topological point of view 

An analysis of the emergent AF helical edge modes 
in terms of topology can be made, but it is less rigor¬ 
ous mathematically than the topological superconductor 
order in the canted AF phase. As shown in Fig. [ljx, 
there is a finite volume in parameter space for which the 
SC contacts are helical metals. The relevant symmetry 
class within the ten-fold way [521 for the infinite system is 
the 2D class D, and its topological invariant is Z, which 
corresponds to the number of chiral edge states at a sur¬ 
face [63] . Within this language the helical SC contact 
has a trivial (zero) Z invariant, since the number of right 
minus left propagating modes is zero. This is actually 
the reason why crossing the 0 = 0[ ns destroys the helical 
edge states without an intervening bulk-gap inversion. 
Therefore, the reason for the existence of helical states 
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for 6 > 0[ ns cannot be found in the standard homotopy 
classification. It is rather an instance of non-trivial valley 
Chern number. 

If one computes the Z invariant in 2D of our sys¬ 
tem, both in the graphene side and on the supercon¬ 
ductor side, one needs to integrate the Berry curvature 
of the Nambu bands. For the superconductor one ob¬ 
tains negligible Berry curvature for all momenta. How¬ 
ever on the graphene side (and choosing the magnetic 
unit cell to compute the bands), one finds that while 
the integrated curvature is zero (hence Z is zero), it 
is the sum of two integer and opposite contributions 
from different valleys. For a specific spin sector (e.g. 
\e t), |h I)), one valley has partial integral 1 and the 
other -1. Therefore, assuming valley symmetry is pre¬ 
served at the graphene/SC interface (that is w = 0 in 
the phase diagram of Fig. [lj^, i.e. the contact is trans¬ 
parent), one can invoke a bulk-boundary correspondence 
principle for each of the two valleys and spin sectors inde¬ 
pendently, which yields two pairs of counterpropagating 
(helical) states (one per valley and spin). If the contact 
is not transparent (w > 0) but intervalley scattering is 
below the threshold w < W[ ns , the helical states will split, 
but the contact will still be metallic (since the splitting 
is smaller than the energy where they cross). Hence the 
AF helical states are an instance of weak topology from 
the two valleys, just like e.g. the helical modes in works 
like Refs. |[4j[65|. Note, however, that the mathematical 


standing of these arguments are less sound than the con¬ 
ventional ones from full-Brillouin-zone invariants, since 
to our knowledge there is no rigorous theorem that guar¬ 
antees the existence of surface states from partially inte¬ 
grated (valley) Chern numbers. 


Appendix E: Square-lattice superconductor 

In Appendix [D] we have considered the superconduc¬ 
tor arising from electrons hopping in a honeycomb lat¬ 
tice and subjected to pairing potential. Nevertheless, the 
fact that the interface states persist even upon doping of 
the superconductor, suggests that their existence goes far 
beyond what our analytic argument might suggest. Ac¬ 
tually, we here show that a honeycomb superconducting 
lattice is not mandatory, so that even an interface be¬ 
tween a square superconductor will give rise to localized 
Majorana states. 

To illustrate this, we show in Fig. [8] the bandstructure 
of an interface between the canted Quantum Hall anti- 
ferromagnet, and the local density of states for a finite 
system. In the case of fully collinear antiferromagnetism, 
the interface sustains a gapless channel. When the mo¬ 
ments are canted, a topological gap in the inter-facial 
bands opens up and localized Majorana modes show up. 
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